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AB3TRACT 

Current  methods  of  shock  design  and  analysis  rely  primarily 
on  the  input  of  large  amounts  of  experimental  data,  assembled  from 
underwater  explosion  tests,  to  describe  the  shock  motion  inputs  to 
their  systems.  The  object  of  this  investigation  is  to  arrive  at  a 
more  complete  characterisation  of  the  shock  motion  input,  a  net 
effect  representing  the  interaction  between  the  water  applied 
pressure  loading,  the  hall  structural  response  and  the  equipment  - 
foundation  response. 

The  model  chosen  is  the  cross  section  of  a  simple 
cylindrical  shell  with  an  inboard  mounted  simple  oscillator, 
immersed  in  an  acoustic  medium.  The  input  to  the  water-shell 
system  is  an  exponentially  decaying  plane  step  wave  incident 
normal  to  the  cylinder  longitudinal  axis. 

The  analysis  employs  the  use  of  four-pole  mechanical 
impedance  techniques  and  Laplace  transform  methods.  An  analytic 
expression  is  derived,  which  has  the  form  of  a  mechanical  quad- 
ripole impedance.  It  relates  shock  wave  pressure  to  equipment 
response,  as  a  function  of  shell  geometry  and  oscillator  mass  and 
natural  frequency.  The  lack  of  rigor  associated  with  the  description 
of  a  distributed  mass  system  as  a  point  phenomenon  along  with  other 
restrictive  conditions,  cause  the  expression  to  be  better  termed  a 
"quasi- impedance1' . 

The  problem  formulation  and  solution  has  limitations  on 
both  the  mathematical  and  physical  side  and  at  this  stage  the  re- 
sults are  to  be  regarded  mainly  as  exploratory.  The  " quasi- impedance" 
is  now  capable  of  predicting  qualitatively  the  effects  of  parameter 
variation  on  the  response  of  simple  equipments  to  shock  wave  inputs. 
Additional  refinements  of  the  shell  model  and  comparison  with  experimental 
data  will  enable  a  more  rigorous  evaluation  of  quantitative  effects. 

Thesis  Supervisor:  John  R.  Baylis 

Title:  Associate  Professor  of  Naval  Engineering 
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A  ■  cross- sectional  area  of  cylindrical  shell 

a  m    radius  of  the  cylinder 

a.^  b,  (i=l>  n)  m    coefficients  in  system  characteristic  equation 

B  ■  non-dimensional  shock  wave  decay  constant 

c  -  sound  speed  in  medium 

e. .  =  general  notation  for  four-pole  parameters 

E  m    Young's  modulus  of  elasticity  for  shell 

F.  •  point  force  per  unit  longitudinal  length 

where  i  refers  to  a  particular  component 

output  or  input  point 
F.  m    modal  component  of  F. 

F!  m     force  per  unit  shell  surface  area 

g.  =  Haywood's  "after  flow"  constant,  average  value 

G(s)  •  general  form  of  transfer  function 

h  •  shell  thickness 
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I  ,  I         ■  incomplete  and  complete  modified  Bessel  fane t ions 

of  the  first  kind  of  order  n,  respectively 

I'  •  I'        ■  first  derivative  of  the  function  with  respect 

no   n 

to  the  argument 
k  p  equivalent  foundation  stiffness  per  unit 

longitudinal  length 
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Le-p  ■  ratio  of  the  longitudinal  length  oi'  shell 

feeling  the  influence  of  the  equipment  reaction 
to  the  actual  equipment  length 

m  m    mass  of  shell  per  unit  area 

ai  -    modal  or  generalized  masses  of  shell  per  unit 

area 

VL  m    mass  of  equipment  model 

n  -    mode  number  and  equivalently  the  number  of 

circumferential  waves  of  mode 
P.  m     free- field  shock  wave  pressure 

I  ■  radiated  and  diffracted  wave  fluid  pressure 

r 

Pin'  *rn       =  G^^nHlzd  or  modal  components  of  P.,  P 

P  «=  shock  wave  peak  pressure 
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a  ,   h  *    generalized  or  modal  deflection  in  the  radial 

direction  and  its  derivative  with  respect  to 
time,  the  generalized  velocity 

6  =  Laplace  transform  of  fi 

r,  9  a.  cylindrical  coordinates 

s  =  complex  frequency  operator 

s.  (1=1,   n)     m    roots  of  characteristic  equation 

t  ■  time 

T  ■  non-c'lmensional  time  -  ct/a 

u  (t)  m  step  function 

U.  =  free- field  shock  wave  particle  velocity 

U  =  radiated  and  diffracted  wave  fluid  particle 
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Uin*  Urn  "    inodal  components  of  U±9  u 

vi  *  point  velocity  where  i  refers  to  the  particular 

component  output  or  input  point 
V^jj  *  modal  component  of  V. 

Vo  -  ?eak  velocity  of  uuaai-terpedance  function 

w  =  radially  inward  90%  velocity  of  shell  as  a 

function  ox 
X(s),  x(t)      m    relative  velocity  between  equipment  mass  and 

shell,  in  complex  frequency  and  time  domain, 

respectively 
zi  ■  general  notation  for  impedance  function 

0  *  shock  wave  decay  constant,  sec 

A  m  notation  for  characteristic  equation 

r  *  gamma  function 

*£  *=  Neumann  Factor;  1  for  n  =  0;  2  for  n  >  0 
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P  *  mass  density  of  medium 
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«-n  ■  shell  stiffness  in-vacuo 


-ix- 


./.'.    '■'-■■  ■ 


'■%'"■'■ 


ito&oift  fifianMX 


V    ■     .      '      '.  :    .    .   J.         .    ■ 


„ 


• 


■ 

Q 


inrf 


potential  function  of  diffracted  and  radiated 
wave3  in  fluid 


General 


function       =  iaplace  transform  of  tne  function  wita  reapect 

to  time 


-\- 


-1- 


I.   IHTBODUCTIGM 

-  ■  -  "i'   ■ 


A.  ahock-Reaistant  Equipment  Design 

In  recent  years  routine  testa  nave  been  conducted  !>y  the 
Navy  to  measure  uhe  resistance  of  naval  snips  to  underwater  shock 
waves  from  non-contact  underwater  explosions-  Th.  ,    rt3  nave  re- 
vealed that  incapacitating  damage  to  shipftcerd  equipments  has  occurred 

at  shock  intensities  of  only  a  fraction  of  those  required  to  produce 

(27) 
lethal  hull  damage*   .  At  the  present  t.me,  considerable  attention 

is  being  devoted  by  various  design  analysis  methods  to  the  prediction 
of  the  responses  that  can  be  expected  by  equipments  to  typical  motion 
inputs.  Because  of  the  complicated  interaction  problems  in  the 
general  case  of  a  shock  wave-shell-equipment  interaction,  shock  motion 
inputs  have  been  arrived  at  only  after  evaluation  of  vast  amounts  of 
experimental  data  from  realistic  underwater  explosion  tests.  Rational 
utilization  of  experimental  data  requires  correlation  between  test 
charge  intensity,  shock  motion  recorded,  and  shell-equipment  inter- 
action parameters,  to  include  relative  mass,  on-board  location,  and 
natural  frequencies ,  By   realistic  modeling  of  equipments  for  which 
dynamic  response  data  from  underwater  exp.lc3ion  tests  are  available, 
it  is  possible  to  solve  for  the  shock  motion  input  at  some  refereno 
point,  i.e.,  inboard  side  of  shell.  A  wide  range  of  dynamic  situations 
have  to  be  investigated  before  a  general  set  of  input  curves  can  be 
arrived  at. 

An  alternate  approach  is  to  tackle  the  shell- shock  wave 
interaction,  making  sufficient  assumptions  to  make  the  problem 
solvable  and  to  investigate: 
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(1)  the  qualitative  validity  of  the  solution 

(2)  possibility  of  the  introduction  of  experimental 
information  to  adjust  the  cuantibttive  solution, 

33ais  thesis  shall  invt:3ti;xte  the  feasibility  of  analytically 
dt 'scribing  the  shock  wars  -  hull  -  equipment  interaction  for  the 
case  of  a  submarine- type  model.  2his  description  will  take  the 
form  of  an  impedance,  Dae  to  non- rigorous,  yet  physically  meaningful, 
assumptions  j  the  result  wiH  be  termed  a  quasi-impedance . 

B.  Shock  Design  Inputs  to  Several  Analysis  Methods. 

A  design  input  may  be  defined  as  a  description  of  the 
shock  forcing  or  motion  function  to  the  design  system.  Dependent  on 
the  system  chosen,  it  could  be  a  velocity  or  acceleration  as  a 
function  of  time  at  some  reference  point  in  the  ship  or  some  pressure- 
time  description  of  the  shock-wave  itself.  A  summary  of  the  design 
inputs  to  the  three  beat-known  design  methods  are  outlined  below. 
1.  Static  Shock  Design  Jfumber 

The  most  rudliiientary  of  the  methods,  the  input  is  a 

shock  design  number,  dependent  upon  equipment  sise,  ship  type  and 

(h-2) 
the  attack  direction*   .  It  ia  applied  to  the  equipment  weight  to 

produce  an  equivalent  force  which  the  foundation  and  mounting  should 

withstand.  A  proposed  modification  to  the      ell  wo -Id  modify 

the  number  by  a  resiliency  factor  based  upon  the  ratio  of  equipment 

natural  frequency  to  the  dominant  frequency  of  the  input  forcing 

function,  as  determined  from  peaks  in  the  frequency  spectrum  of  the 

input  motion*  '. 
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2.  Normal  Mode. 

The  normal  mode  approach  currently  describes  an  equivalent 
shock  motion  input  (frequency  domain)  to  a  as fchematical  model  at 
some  inboard  fixed-base  reference  point,  $.«.,  for  I  submarine,  the 
inboard  side  of  the  pressure  hull.  The  fundamental  description  of  the 
input  is  called,  the  shock  spectrum.  The  shock  spectrum  is  defined 
as  the  trace  of  the  maxima  absolute  response  of  a  single  degree 
of  freedom  system  to  the  applied  rigid  base  nation  as  the  freque/. 
of  the  qjatlllt  a*"  If    -.   from  0  -  *•  Its  gptntiXJi  as  a  design 
input  for  naval  shock  lies  in  the  fact  that  it  fails  to  account  for 
structural  interaction  between  the  oscillator  and  the  foundation 
exciting  it,  The  spectrum  provides  no  information  with  regard  to 
time  history,  but  rather  an  upper  bound,  on  the  maximum  value  of  the 
response  as  a  function  of  frequency.  sTwtimlim  of  uhe  non-damped 
vibration  absorber  from  classical  vibrations  *  '  and  discrete-mass 
laboratory  experiments  have  led  some  oh  •  .      to  propot, 

that  the  values  of  interest  in  a  shoex  specurum  tend  to  lie  in  MM 
©ills  of  the  plot,  rather  than  at  'the  peaks.  A  fiducial  en/ elope  of 
the  peaks  is  valid  as  an  input  when  dynamic  reaction  of  the  driven 
component  is  unable  to  affect  the  motion  of  the  driving  component. 
So  as  to  account  for  the  interaction.,  the  fundamental  she,     sta 
is  modified  by  equipment  weight,  on-board  location,  and  natural. 

frequency.  The  result  is  termed  a  "design  shock  spectrum"  smd  is 

(.2 

the  design  input v     ' .     By  assembling  enough  ital  d~ 

(equipment  response)  ranging  over  a  sufficient  banc!  of  equipment 
natural  frequencies,  relative  mass  ratios  between  equipment-foundation 
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aad  bull  mass,  on-board  location,  charge  (shock-wave)  intensities, 
and  ship  types,  design  inputs  can  be  as-rived  at.  Each  point  on  a 
design  input  chart  requires  equipi^ut-founG-^ion  modeling  and  inverse 
normal  mode  solution  for  the  forcir-g  function  required  to  produce 
the  measured  response.  For  the  submarine  case,  a  modification  to 
this  approach  has  been  suggested  which  moves  the  inference  ooint 
outside  the  shell  and  adds  hull  flexibility  to  foundation  flexibility, 
Hull  resiliency  would  be  experimentally  determined*""  M   . 

3.  Iteration. 

Essentially  a  numerical  method,  it  describes  the  hull  • 
foundation  -  equipment  interaction  by  a  non-linear  two-mass  system. 
It  is  solved  by  convergent  iteration  on  the  equations  of  motion. 
The  design  input  is  a  pressure-time  history  of  the  free -water  shock 
wave  which  is  the  forcing  lunetion  on  a  discrete  mass  describing  the 
hull  "effective  mas3"w  ;. 

C.  Mechanical  Impedance  Approach. 

Concurrent  with  renewed  interest  in  the  sdaptation  of 
classical  dynamics  to  complicated  structures,  has  been  the  adoption 
of  an  analytic  tool  from  electrical  engine ering  known  as  the 
impedance  technique. 

The  use  cf  impedance  permits  combinations  of  a  large  com- 
plicated arrangement  cf  components  to  be  handled  as  a  single  entity. 
Computational  efficiency  is  increased  "because  only  input  and  c 
data  are  handled,  eliminating  explicit  introduction  of  all  the 
parameters  characterizing  each  component  "black  box" .  Further,  the 
parameters  contained  in  each  "black  box"  are  a  function  only  of  th 
particular  component  and  are  not  influenced  "by  the  component* 8 
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surroundings.  The  use  of  such  a  technique  la  the  analysis  of  shock 
motions  appears  Ideal  in  that  it  provides: 

(1)  a  flexible  procedure  allowing  for  the  introduction 
of  empirical  and  quasi-analytic  functions  describing 
"black  box"  interactions. 

(2)  direct  utilization  of  experimental  data 

(3)  maximum  ease  of  evaluation  of  effects  on  the  system 
response  due  to  component  changes,  i.e.,  removal, 
alteration,  replacement. 

To  characterize  the  interaction  of  a  "black  box"  with  its 
surroundings,  the  impedance  technique  uses  term3  such  as  load  or 
driving  point  Impedance,  output  impedance,  and  transfer  impedance. 
A  further  refinement  of  the  impedance  technique  formalizes  these 
interaction  impedances  into  two  linear  equations  which  relate  the 
input  point,  1,  to  the  output  point,  2,  in  terms  of  their  four-pole 
parameters: 

Fl  =  ell  F2  +  e12  V2 

Vl  "  e21  F2  #  e22  V2 
This  approach  is  called  the  Four-Pole  Parameter  Technique  and  is  an 
analytic  tool  in  this  investigation^  *™*, 

D.  Shock  Wave  -  Bull  Interaction 

Pressure  wave  interaction  with  cylindrical  obstacles  was 
investigated  classically  by  Raleigh*   .  The  problem  has 
received  revived  interest  in  the  application  of  the  basic  approach 
to  obstacles  whose  thin  cross  sections  are  allowed  to  deform 
elastic ally. 
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Tlie  earliest  of  these  studies s  ' '  -used  a  modal  analysis  and 
retained  the  full,  exact  form  of  the  diffracted  and  radiated  wave. 
The  analysis  was  handled  by  transform  methods.  Due  to  the  complexity 
of  the  solution  (presence  of  complex  poles  and  branch  points),  inversion 
by  contour  integration  for  more  than  several  modes  was  tedious.  Mindlin 

and  Bleicir  ',  by  introducing  an  approximation  valid  for  small  time, 

(V7) 
simplified  the  solution  of    and  presented  numerical  results  for 

the  first  three  modes  of  cylinder  response.  The  simplification 

amounted  to  giving  the  propagation  characteristics cf  &  plane  wave  to 

the  radiated  and  diffracted  wave.  Extenslonal  effects  in  all  modes, 

(15) 
but  the  dilatational  mode,  were  neglected.  Barcn  *■'   extended  the 

analysis  of    'to  include  both  inextensional  and  extensional  effects 
in  all  modes,  using  the  plane  wave  approximation.  More  recently, 
Haywood   '  has  contributed  a  derivation  allowing  for  a  more  accurate 
approximation  to  the  characteristic  impedance  of  a  cylindrical  wave 
allowing  for  an  "after- flow*  effect.  The  approximation  reduces  to 
that  of*  '  for  the  case  where  the  modal  "after- flow''  constant 
reduces  to  zero.  This  treatment  extends  the  time  validity  range  to 
times  of  the  order  of  shock  wave  transit  time,  while  reducing  accuracy 
for  small  values  of  time. 

These  investigations  have  been  restricted  to  the  two- 
dimensional  problem,  i.e.,  response  to  long  shocx  waves.  Longitudinal 
effects  in  the  absence  of  cross  sectional  deformation  are  investigated 
in^33)  aaflW).  For  ^iQ   'whipping"  mode,  the  model  is  a  cylinder  with 
rigid  cross  section,  which  is  allowed  to  bend  eiastically  about  its 
longitudinal  axis  in  response  to  a  local  (spherical)  step  wave.  For 
the  case  of  a  symmetrical  cylinder-wave  situation  (plane  wave  -  cylinder), 
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■  • 

. 

-  ■■.da>0'   9V4U1    fcO. 

«rt     .»y*v  q»*t  (XaoliMfcp)  i*oa£  jb  o*  »&aoqw«  cl  ,X 


-7- 


the  results  of  Murray w  '   compare  in  general  with  the  tranalatioaal 
(n  =  1)  mode  of  MindLin  -  Bleieh^l8\ 

£.  la- Vacuo  Shell  Kormal,  Modes* 

The  natural  frequencies  for  a  circular  cross- section 

in-vaeuo,  as  applied  in  each  of  the  previous  approaches,  are 

(2) 
derived  utilising  LaGrange's  equations  by  T5jaoshenkox  ',  The 

in-vacuo  frequencies  of  vibration  of  a  simply  supported  thin  cylindrical 

shell,  reinforced  by  equally-spaced,  circular  ring  stiffeners,  are 

derived  lay   Galletiyv  J  .  Bleictr  **   has  derived  an  approximate 

method  for  calculating  the  normal  modes  of  ring- stiffened  cylinders. 

found  to  agree  with  the  more  exact  treatment  of   '  within  4  10  -per 

cent.  These  methods  satisfactorily  describe  the  natural  frequencies 

of  a  more  realistic  model  for  lower  modes  without  undue  labor, 

(45) 
The  general  form  of  the  results  of    '  are  listed  in  Appendix  D. 
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II.  THEOBETICAL  PROCEJXIRE 

A. Description  of  System  Model. 

The  hull  is  characterized  by  a  simple  cylindrical  shell 
with  a  circular  cross*  section  immersed  in  an  infinite  acoustic 
fluid.  The  response  of  the  shell  will  be  elastic.  Mounted  at 
Point  A  on  the  shell  is  a  simple  oscillator  representative  of  a 
foundation  -  equipment  chain,  whose  elastic  response  is  to  be 
determined.  The  equipment  response  is  defined  to  be  one- dimensional, 
Although  the  physical  situation  is  defined  to  be  two-dimensional 
(Fig.  1),  longitudinal  effects  are  compensated  for  in  part  by  the 
use  of  an  empirical  factor,  called  the  effective  length  ratio. 
No  "whipping"  effects  are  considered.  The  equipment  -  foundation 
response  time,  while  remaining  a  variable  in  the  solution  quasi- 
jbtapedaace,  is  stated  to  be  of  the  order  of  the  shock  wave  transit 
time  (5-6  msec.).  The  shell-simple  oscillator  combination  is 
subjected  to  an  acoustic,  plane t   exponentially- decaying  step  wave, 
normally  incident  to  the  longitudinal  axis  of  the  cylinder.  The 
use  of  linearized  theory  i3  valid  only  for  pressure  waves  which 
are  not  too  intense  and  boundary  excursions  which  are  relatively 
small.  Symmetric  modal  deformation  of  the  cylindrical  shell  is 
assumed. 

B.  General  Approach. 

The  sequence  of  events  in  the  derivation  of  the  cylinder 
quasi- impedance  will  be  outlined  before  starting  the  mathematical 
details . 


-- 
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1.  Define  the  shell  response,  the  system  input  force 
and  velocity,  arid  a  potential  function  for  the 
radiated  wave  in  terms  of  Fourier  expansions  of  their 
generalised  coordinates,  generalized  force  and 
velocity,  ana  modal  potential  function. 

2.  Define  an  effective  modal  raass  of  the  shell. 

3.  Setup  of  the  elation  of  motion  of  a  point  on  the 
shell.  The  shell  is  immersed  in  the  medium  and  the 
point  feels  the  interaction  with  an  equipment  - 
foundation  chain  mounted  inside ,  The  shell  is  act 
upon  "by  the  summation  of  the  incident,  radial 

and  diffracted  pressure  waves. 
k.     Description  of  the  natural  modal  frequencies  and. 
shapes  in- vacuo  for  the  shell  model  chosen. 

5.  Approximate  the  characteristic  impedance  of  the 
radiated  wave. 

6.  Prescribe  necessary  boundary  conditions. 

|.  Using  (5)  and  (6),  eliminate  the  fluid  potential 

function  from  the  modal  equation  of  motion. 
Q.  Express  the  modal  equation  of  motion  in  terms  of  :, 

Laplace  trausfonas  and  convert  to  normal  Tour-pole 

parameter  format. 
9.  Set  up  the  four-pole  parameters  for  an  appropri; 

model  of  the  .foundation  -  ecj.apment  chain. 
10.  Using  matrix  techniques.,  combine  the  shell  and  the 

equipment  -  foundation  jbnpeda&ec  fUQctiCBcus 
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11 .  Invert  the  resulting  Laplacian  expression  of  (10) 
to  the  time  domain  solving  for  the  equipment  re- 
sponse to  one  mode  of  the  shell  response  to  the  input 
shock  wave. 

12.  Sum  the  responses  of  the  equipment  to  the  modal 
excitations  of  the  shell  for  as  many  modes  as 
required. 

13.  Convert  to  response  relative  to  hull  if  desired, 

C«  Impedance  ffirametera  . 

A  general  discussion  of  the  mechanics  of  forming  cuadripole 

(23) 
parameters  for  a  wide  variety  of  physical  systems  is  available  in    . 

The  brief  discussion  included  here  is  solely  for  the  sake  of  continuity, 

1 .  Definitions 

A  linear  system  which  possesses  a  single  input  point.  I, 
and  a  single  output  point,  2,  may  be  described  by  a  pair  of  linear 
e<?uations  of  the  form 

h  *   eU  F2  (  e12  V2  W 

Vl  "  e21  F2  *  *22  V2 
where  e...,  e.p,  e^,  and  e^p  are  the  quadripole  or  four-pole 

parameters  of  the  system.  Solving  Equation  (l)  for  the  four  poles 
gives: 

e1,  =  — -  with  V,-  ■  0,  i.e.,  station  2  blocked; 
e12  "  V~  ^^  F?  =  Q'   i-6-.  s'cation  2  free: 
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e01  =  rr  with  V0  =  0,  i.e.,  station  2  blocked; 

21   F2       2    ' 

e22  =  V^  vl^il  ^o  "  °>  i°e.,  station  2  free. 

The  e^,  parameter  is  thus  the  force  transmissibility  from  station  1 
to  station  2  where  station  2  is  "blocked.  Parameter  e,2  is  the 
transfer  impedance  between  station  1  and  station  2  with  station  2 
free-  Parameter  e^,  is  determined  by  the  transfer  impedance 
between  stations  1  and  2  with  station  2  blocked.  Parameter  eg2  is 
defined  in  terms  of  the  velocity  transmittance  from  station  1  to 
station  2  with  station  2  free.  Equations  (1)  are  conveniently 
expressed  in  matrix  notation  as: 

[?]  -  rl11 112']  I"  v2i  w 

LV   l-e21  e22J  LV 
The  first  bracketed  term  on  the  right  is  called,  the  quadripole 
parameter  matrix.  Matrix  algebra  provides  a  formalized  and 
practical  method  of  forming  component  quadripoles  from  fundamental 
elements,  and  of  obtaining  the  overall  system  quadripole  from  its 
components . 

2.  Qmdripole  Impedances. 

The  elements  of  the  four-pole  matrix  are  naturally 
related  to  the  various  types  of  mechanical  impedance  -hat  describe 
the  component.  The  relationship  between  the  two  are  listed  below 
in  order  to  provide  correlation  between  the  two  methods  of 
nomenclature . 

The  load  or  driving  point  impedance  characterizes  the 
elements  that  follow  it: 


. 


■ 


- 


■ 

- 


. 


.U-J 


i»£rxo  oi 


M 


> 


-13- 


The  input  or  driving  point  impedance,  defined  by 

K 

*  vi 

is  dependent  on  both  the  quadripole  and  on  its  load  Impedance.  It 
is  related  to  the  load  and  the  four-pole  "by 

z  ~  e-u  z3 *  eig 

1   e^  Zg  .  e22 

The  input  transfer  impedance  is 
h 

ha  "  %  -  ~n  22  *  £i2  * 

The  force  transfer  function  is 


Fx 


e-jp/Z,. 


The  velocity  transfer  function  is 

VV2  *  e22   e21  h 

3.  Forming  the  Parameters. 

For  some  simple  models  direct  combination  of  the  known 
parameters  for  fundamental  elements  (Appendix  C)  using  matrix 
algebra  will  produce  the  quadripole  for  the  system.  The  definitions 
given  above  also  provide  a  means  of  introducing  experimental  de.\ 
into  the  analytic  solution.  However,  in  che  general  case  it  is 
necessary  to: 

(1)  define  the  system 

(2)  describe  the  system  in  tens  of  a  single  pair 
of  input  and  output  junctions 


zl 


fix' 
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(3)  odve  -tiie  performance  equation  using  operational 

calculus  methoda 
(k)     resolve  the  equations  into  the  form  of  equation  (1) 

or  (2). 

D.  Analysis  or  Shell* Shock  Wave  -  Foundation  -  Equipment  Interaction 
1.  Incident  Wave. 

The  input  is  defined  as  an  exponentially  decaying  step 
wave  of  the  form 

fc.  -  t  e"3t  u(t). 
i   o 

The  simple  step  is  translated  across  the  cylinder  with  the  use  of 

a  time  variant  delay  factor 

t  »  -  (  1  -  cos  9)    (origin  taken  at  cylinder  axis) 

C 

so  that 

Pi  -  PQ  ■***  '  T)  u  (t  -  t)  (3) 

The  total  pressure  and  the  radial  component  of  the  particle  velocity 
are  expanded  into  Fourier  series  of  the  form 

Pi  ■  I  pin  cos  a 
n=Q 


Ui  =  I     Uin  cos  n  ^ 

n-0 

Tne  Fourier  modal  coefficient  P.  is  obtained  from  equation  (h)  by 

modal  orthogonality  after  multiplying  right  and  left  sides  of 

equation  {k)  by  cos  n  I,  integrating  iron  0  to  gffj  interchanging 

the  order  of  3uiaaatian  and  integration.  The  integrand  is  zero 

from   to  Jfe  by  the  definition  of  the  step,  u(t  -  t). 
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This  procedure  yields 

it) 

Pin  =  n  J  Pi  (  >t}  co3  r'    °  <  r'   (t)  <  * 
0 

•where   (t)  is  the  angular  position  of  the  incident  wave,  Using 

the  Laplace  time  -  lag  theorem  for  translation  in  the  tiaae  domain -> 

substituting  from  equation  {k)   and  transforming  with  respect  to  the 

non-dimensional  time  T  =  ct/a  yields, 

-8        cos"1  (1-T) 

0  (5) 

Similar  development  applies  to  the  radial  modal  particle  velocity. 

Introducing  an  additional  cos  I  because  of  its  radial  definition 

and  using  the  characteristic  impedance  of  a  plane  wave,. 

(t) 

Uin  *  -  *  ^C     Pi  '*'*'  COS  n  °  k©9   d    0  <^  -^  2rr 
0 

and 

-s   (t> 
(s  i-  B)©^  «  i  x  ^  JV  8  cos  '  cos  a  l)  cos  |  d   0  ^   * 

0  (6) 

At  I  =  2a,  the  integral  in  equation  (5)  becomes  a  definition  of 
the  modified  Bessel  function .  During  transit  time,  they  are 
referred  to  as  incomplete  modified  Bessel  functions.  Their  numerical 

j.uaticn  is  discussed  in  Appendix  B*  Using  this  definition, 
equations  (5)  and  (6)  reduce  to: 
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-8       P 

Jin   3^B  n  c  x  n  <J) 

(prime  denotes  derivative  with 
reapect  to  s.) 

2.  Interaction  Pressure  Field. 

Contact  between  the  incident  wave  and  the  cylinder  pro- 
duces scattered  and  radiated  acoustic  waves  whose  potential 
function  must  satisfy  the  wave  equation.  The  scattered  wave 
arises  because  of  the  presence  of  the  cylinder.  2$te  radiated  wave 
Ml  caused  by  the  response  motion  of  the  cylinder .  The   scatter 
and  radiated  waves  are  assumed  to  move  radially  oat  from  the 
cylinder.  She  radiated  pressure  and  velocity  are  defined  in  terms 
of  the  potential  function  I  by 
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r     r 
Modal  contributions,  corresponding  to  the  in-vacuo 
generalized  coordinates  of  cylinder  normal  modes,  are  obtained 
by  expanding  the  velocity  potential,  thus  P  and  U  into  a  Fourier 
cosine  series  in  !  of  the  form  of  equation  (4).  Equating  coefficients 
of  the  expansions;  yields: 
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where  the  generalised  potential  function  also  satisfies  the  wave 
equation  in  cylindrical  coordinates 

n   1  %   _      1     n 

ar*»F'i    n  r: 7  -rr 

<fP  r        C    :^t 

3-  Cylindrical  Wave  Approximation. 

The  use  of  the  complete  form  of  the  wave  equation  above 
to  provide  a  relation  between  P   and  U   (characteristic  impedance) 
yields  a  solution  whose  inversion  from  the  complex  frequency  domain 
is  prohibitively  tedious. 

An  approximation  by  Haywood*  '  to  the  characteristic 
Impedance  of  a  cylindrical  diverging  wave  yields  results  valid 
for  times  of  the  order  of  the  shoclc  wave  tranfiit  time,     Haywood's 
derivation  yields: 

(8) 


or 


r^r    '  c>5  "  r  n 


*u    ,  if    g„ 
_rn   1^   rn   to  p 

it  ~  pc  51   pr  ra 
The  modal  constant  eL  is  the  average   ft  .ion  whose  limits ,   while 
time  ranges  between  0  and  *>,  is: 
0  <  g 

1/2  <  gQ  <  n     f or  n  >  1  {-)) 

The  main  effect  of  the  approximation  i~  i  «  fln:~     -**%', 

or  "afterflov "  to  the  fluid,  following  the  passage  of  the  pressure 
wave,,  which  goes  to  zero  only  after  the  rressur .  hau  fallen  below 
its  equilibrium  value.  For  distances  infln.     'ar  from  the 
cylinder,  eviction  (B)  reduces  to  the  plane  wave  appi-cximation. 


(9»- 


IVSW   I*oJtX 


arf*  »di 


o«ft«  (ft)  «oj  ,/x9Aa 


-18- 

k.     Natural  frequencies  of  Shell  Model. 

For  s         tee,  the  simple  circular  unstiffeaed 
shell  model  is  considered  adequate.  Extensioa&l  effects  in  all 
modes  except  the  dilatations!  (a  *  0)  are  neglected.  The 
derivation v  ' v  'of  the  normal  mode  frequencies  for  a  simple 
circular  cross  section  yields: 

E 


with  E  - 


(1  -„2) 


■a-  2  (taexteasioaal  -  ilexural)  ■  n  ^  "  a  i  x  -  x  |  x  h 
a  1  h  n^     ps   A   aA 


The  natural  frequency  of  the  translations!  mode  a  =  1  is,  of  course, 
equal  to  zero. 

In  order  to  reduce  the  order  or  the  characteristic 
equation  for  a  later  approximation^  we  now  define 

(«  •  -=  as  the  modal  shell  stiffness, 
a   hl 


where  la  *  m,  (n  ■  0)   and 


n 

>.  Equation  of  Motion  at  Point  on  Shell. 

Summation  of  the  modal  forces  acting  at  point  A  in 
Fig#  1  gives: 

at" 
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where  (L  ia  s  generalized  coordinate  and  is  defined  in  terms  of 
the  total  radial  deflection  w  by: 

W  "  )    %i    COS  a  6 

6.  Boundary  Condition. 

At  any  point  on  the  surface  of  the  shell,  the  modal  velocity 
of  the  shell  must  be  equal  to  the  net  particle  velocity  of  the  in- 
cident and  radiated  -  diffracted  modal  waves 

*n  -  uin  "  um  (l3) 

7.  Ret  Shell  Response  Function. 

An  operational  flow  chart  combining  the  preceding  steps 
is  given  in  Figure  2. 

Transforming  equations  (8),  (12)  and  (13)  by  Laplace 
transforms  with  respect  to  time  yields 
2   -2 

<VV  "  ">  I  ■  »!.  -  "a*  ^  <»!.  -  V    <1U> 

a 

Rearranging  and  changing  the  time  variable  to  a  non-dimensional 
time  T  (=  ct/a)  yields: 

+  SS  pC  UjJ  (15) 

where  A,  the  characteristic  equation  of  the  shell  response,  equals: 

C  PC 
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Fram  equation  (7), 

f*»     -  v> (s) 

therefore,  equation  (i>)  reduces  to: 

Tl 

pc  (A)  ^  =  -s  (s  +  g^  F^  ♦  J  s(s  +  ijj)  +  s2  y^I  p. 


la 

n 


(1/) 
Introducing  the  definition  of  P.  for  step  decaying  shock  wave 
equation  (7),  we  obtain 


T> 


I  (s  +  B)(s  ♦  Zq) 


«"s  L- 1° +  ^)  v  <s>  * *a  *i  <•>. 


F. 


in 


, pO  A  (o  i  B) /lAv 

L3  (s  H  sn)  **  (s)  *  8<"  *n;    (fl) > 

E»  Equipment  -  Foundation  and  Shell Four-Pole  Parameters. 

In  general;  there  is  no  limit  to  the  complexity  of  the 
model  of  the  equipment  -  foundation  chain  mounted  at  A  in  Figure  1, 
except  that  it  be  capable  of  representation  by  point  attachment. 
For  the  purpose  of  this  investigation,  a  simple  oscillator  was 
chosen   The  mass  and  natural  frequency  zaay  be  varied  to  d«  treriulne 
their  effects  on  the  shell  -  equipment  interaction* 

The  mass  four-pole  as  derived  in  Appendix  C  is: 

1   stt,  1 
_0   1  *] 

The  spring  four-pole  is  found  to  be; 

f'l    CT 
Ls/k   1. 

They  are  connected  in  cascade,  thus  aj-e  combined  oy  matrix 

multiplication  to  give: 
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This  four-pole  has  been  transformed  with  respect  to  time.  In 

order  to  conform  to  equation  (13),  we  change  the  time  variable 

to  non-dimensional  form  T   (=  t-).  The  four-pole  beccanes: 
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ye  see  that  elation  (l8)  fits  the  general  form  of  one 
of  the  linear  four-pole  equations  (l).  There  are  several  points 
to  be  discussed  or  rationalized  before  we  give  (l8)  an  impedance 
tag. 

(1)  Equation  (15 )  describes  the  modal  response  of  a 
distributed  parameter  system  to  the  spatially 
distributed,  time  variant  force  of  the  pressure 
loading  and  a  dimensionally  similar  equipment  reaction 
force i  Implicit  in  the  impedance  concept  is  the 
necessity  for  point  input  -  output  description  of 

the  system.  The  application  of  the  impedance 
concept  can  be  rationalized  if  we  consider  as  the 
input  to  point  A  en  equivalent  point  force  which 
would  produce  the  same  response  at  A  as  the  distributed 
pressure  loading.  Attached  to  the  output,  of  point  A 
is  the  equipment  model. 

(2)  To  further  satisfy  our  point  requirement  for  the 
cross  section,  we  stave  that: 
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*2n  (ecIuation  (4.5)  )  =  -^  (2Jl) 

vhere  jf'  is  an  equivalent  point  force  for  the 
cross  section,  but  remains  a  distributed  force 
longitudinally  representing  the  equipment 
reaction  force  per  unit  of  longitudinal  length. 

h    "    **  \ff  (22) 

where  L   is  defined  to  be  the  ratio  of  the 
longitudinal  length  of  shell  feeling  the  influence 
of  the  equipment  reaction,  to  the  actual  equipment 
length   L  fif,  is,  of  course,  an  empirical  input. 
It  vlli  be  a  function  of  relative  stiffnesses  and 
masses,  but  a  range  of  2-5  is  considered  reasonable 

(3)  The  interaction  function  %M  defined  here  only  for 
an  exponential 3 y  decaying  step  wave.  However; 
physically,  this  is  quite  acceptable,  since  by 
superposition  and  the  use  of  the  Laplace  delay 
factor,  the  characteristics  of  any  normal  pressure 
loading  con  be  simulated. 

(k)    The  net  impedance  function  is  defined  in  terms  of 
an  infinite  series  of  modal  Impedances. 

(5)  The  complexity  of  x&e  model  prevents  coMplete 
description  of  the  four-pole  parameters*.  The 
problem  is  in  essence  similar  to  the  transmission 
line. problem  in  electrical  engineering*  There ; 
plane  waves  travelling  along  a  transmission  line, 
described  hy   a  characteristic  impedance ,   meet  an 
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interface  or  load  impedance  different  from  that  of 
the  transmission  line.  Due  to  the  interaction, 
standing  waves  are  set  up  in  the  transmission  line 
altering  the  voltage  -  current  relationship  (input 
impedance).  The  analogy  is  made  by  regarding  the 
transmission  line  as  the  medium  surrounding  the 
cylinder ,  and  replacing  voltage  and  current,  by 
force  and  velocity-  The  essential  difference  and 
reason  that  the  analogy  can  be  extended  no  further 
arises  from  the  fact  that  with  the  cylinder  the 
problem  is  no  longer  bilateral.  The  input  wave 
impedance  differs  from  that  of  the  output. 
In  order  to  form  the  four-pole  parameters  according  to 

the  standard  recipe 3   we  would  first  define  the  system  components 

and  component  point  inputs  and  outputs  ass 
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b)  Input  transfer  impedance  for  the  water  component  is: 
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c  Force  transfer  impedance  for  the  water  is: 


d.  Velocity  transfer  impedance  for  the  water  is: 
Ik 
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e,  Load  impedance  at  the  'water  output  function  is: 

'a 

It  is  possible  to  describe 
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using  equations  (lt>),  {?),  (13),  but  we  have  not  adequately  described 
the  interaction  to  be  able  to  define  both  without  obtaining  one 
redundant  set*  However,  for  cases  where  the  response  of  the  ||M 
element  in  a  chain  is  desired,  i.e.,  output  end  is  free,  force  equal 
to  zero,  one  of  the  two  four-pole  equations  13  adequate  to  describe 
the  motion  response  All  interactions  enroute  are  still  accounted 
for.  Such  is  the  case  in  this  example- 

Thus,  we  have  not  described  tl        components  entirely - 
In  fact,  we  have  shunted  the  complicated  Interaction  at  the  shell 
surface  and  have  dealt  with  a  system  whose  limits  are  defined  at 
the  input  ~by   the  free  stream  shade  wave  anc  at  the  output  by  the 
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inboard  junction  to  the  equipment  -  foundation  four-pole. 

With  these  rather  restrictive  conditions  governing  our 
transfer  function?  we  call  it  a  quasi- impedance. 

The  shell  Junction  equation  (i6)  is  now  combined  with 
the  four-pole  of  (20)  and  definition  equation  (21)  to  produce  after 
some  rearrangement; 
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(25) 
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(27)  using  the  following  recurrence  relation: 

*a  («•)  *  ||  *a  («•)  ■  •  Vl  k»}  -  f  Ia  (as)  (2o) 

the  velocity  is  inverted  from  the  complex  frequency  domain.  Use 
of  the  residue  theorem  from  complex  varialale  theory  (4)  makes  the 
procedure  quite  straightforward. 

3!he  velocity  response  for  all  modes  other  than  translations! 
(n  =  l)  takes  the  form: 

V3  (T)  =  ~(s1  +  B)(s2   -  B)(s3  +  B)(%  ,   B)(Sp^r~ 

+  «n  vo  }  ~ — i: ~~ —      " 
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where  B^  •  (^  +  BXa^  -  S2)(3X  -  33^3i  "  s2i^si  "  O 

H2  •  (s2  ♦  B)(s2  -  s1)(sg  -  s3){82  -  \)(*2  *  O 
B^  *  (s3  ♦  B)(s3  -  s1)(s3  -  s2)(s>3  -  %)U5  -  »5) 
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^  =  (»5  *  B){s5  -  3x)(s5  -  s2)(g5  -  s3)(s5  -  ak) 

la  the  traasia clonal  mode  (a  •  1),  the  shell  stiffness., 

-2 

•L*  goes  to  zero.  Tae  coefficients  of  the  characteristic  equation, 

9  (s^  +  a^  s  +  a  9  +  ag), 
reduce  to 

%  '  *a  J' •  § 

a3  "  7  <  W  *  \   %  (30) 

&2  /-  2        -  .  fiWch      1      k      -   , 

a2  =  J|  "Co*  Sa  *  ps  1^-  "iq    ^    V 

aQ»0 

The  calculation  of  the  residues  mist  be  repeated  due  to  the  double 
pole  at  the  origin,  The  translational  velocity  of  the  oscillator 
aass  is  found  to  be  of  the  form: 
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where  E^  •   (a1  -  Sg)^  *  s3^ 

H2  -  (s2  -  3;l)(32  -  s3) 

H3  -  {8X  -  b3)(32  -  a3) 

The  velocity  of  the  equipment  mass  as  derived  here  is  in  an  absolute 
reference  space.  The  absolute  velocity  is  a  satisfactory  measure 
of  the  effects  of  system  and  component  parameter  variation.  However, 
from  the  standpoint  cf  equipment  damage  criterion,  a  more  useful 
quantity  is  the  relative  motion  between  the  hull  and  the  equipment. 
Relative  motion  will,,  of  course,  govern  the  maximum  stress,  deflection, 
and  acceleration  that  the  equipment  and  foundation  will  see.  The 
power  of  the  impedance  technique  Is  revealed  in  the  simplicity  with 
which  the  relative  velocity  is  calculated. 

(1)  Equation  (27)  gives  the  equipment  modal  velocity 
in  frequency  (s)  domain, 

*3n  (a)  "  ofsl 

(2)  Equipment  -  foundation  four-pole  parameters  given  as: 
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(3)  Since  F^  -  0, 
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(U)  Therefore, 
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(5)  Inversion  to  time  domain,  if  desired,  is  identical 
in  procedure  to  that  of  V-^  (a),  njhe  poles  of  the 
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response  function  are  of  course  identical  since 
the  characteristic  equation  for  the  system  is 
lefcango; 
(..)  The  relative  velocity  is  then: 

X   (s)  •  v^  (a)  -  v^  (s)   or, 

x  (T)  -  V^  (T)  -  V^  (T).  (33) 

F.  Numerical  Procedure. 

Equations  (29)  «nd  (31)  will  form  the  basis  of  the  time 
response  calculations.  While  the  entire  procedure  is  straightforward, 
the  number  of  cycles  required  to  trace  the  time  history,  coupled 
with  the  superposition  of  an  indeterminant  number  of  shell  mode..;, 
and  the  large  amounts  of  complex  number  arithmetic,  lead  to  the  use 
of  a  digital  computer.  Since  the  general  form  oi  she  residues  io 
determined,  the  machine  needs  only  to: 

(1)  accurately  determine  the  roo~s,  real  and  complex, 
of  the  characteristic  equation. 

(2)  perform  numerical  quadratui'e  or  series  approximation 
(Appendix  B)  to  the  incomplete  Bessel  functions  at 
steps  sufficiently  small  to  provide  time  history 
desired 

(3)  evaluate  equations  (29)  and  (31)  as  functions  of 
time 

(k)     superpose  modal  responses. 

The  frequency  response  (s  »  ia>)  of  the  quasi-impedanee 
is  readily  available  by  choosing  the  range  sat  discrimination  necessary 
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Once  ■iitiiftM>tc»|^y  prograoaiaed,  ill  paramsterj  can,  of 
course^  be  varied,  using  systematic  procedures  to  determine 
existing  trends.  For  the  initial  computer  rons,  typical  parameters 
were  arbitrarily  chosen  to  characterize:  the  shell  model,  and  shock 
wave.  Three  separate  combinations  of  equipment  -  foundation  model 
mass  and  natural  frequency  were  ehosen  in  a  typical  range  to  order 
to  evaluate  the  ability  of  the  quasi-tai     j  to  describe  interaction 
effects . 

So  as  to  be  able  to  check  by  hand  initial  calculations 
using  a  closed  form  root-finder,  the  order  of  the  characteristic 
equation  was  reduced  from  quintic  to  euartic.  This  reduction  was 
accomplished  by  neglecting  the  shell  mass.  The  maximum  error  will 
be  experienced  for  times  shortly  after  the  shock  wave  arrival. 
The  natural  periods  of  the  equipment  models  chosen ;  are  of  the 
order  of  transit  time.  However,  it  is  obvious  from  equation  (2  ) 
that  the  complete  interaction  is  a  function  of  relative  hull  - 
equipment  masses.  Therefore,  no  physical  justification  is  offered, 
but  rather  the  argument  that  the  general  ability  or  inability  of 
the  function  to  characterize  the  interaction  will  be  retained. 
Details  of  the  programming  procedures  are  listed  in  Appendix  A. 
The  computer  language  is  IBM  FORTRAN  and  the  computer  IBM  7090. 
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Ill,     DISCUSSION  OF  HE3UIff3 
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A  quasi- inpedance  function  has  been  derived  .Cor  •  simple, 
but  physically  reason     W4t1.  The  function  possesses  flexibility 
for  extension  to  analytic  models  of  a  move  complex  nature  without 
changes  in  its  basic  form.  It  further  provides  an  opportunity 
for  incorporation  of  empirical  data  in  the  form  ©f  measured 
component  impedances  and  shell  geometry  variables  into  an  analytic 
description  on  a  rational  "basis.  The  faction  describes  the 
equipment  model  response  in  terms  of  physically  meaningful  parameters: 

a.  equipment  -  foundation  mass,  stiffness,  and  natural 
frequency  or  experimentally  determined  impedances 

b.  shell  mass 

c.  shell  in-vacuo  natural  frequencies 

d.  shell  geometry 

1.  radius 

2.  cross-section  moment  of  inertia 

3.  shell  thickness 

k.     f raise  spacing,  area  and.  moment  of  inertia 

with  the  use  of  Appendix  D 
5.  empirical  input  of  "effective  length  ratio* 
e-  3hock  pMNi  characteristics 
1.  peak  pressure 

exponential  decay  constant. 
Complete  evaluation  can,  of  course  he  made  only  after 
completion  of  a  program  of  systematic  parameter  variation  in 
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conjunetion  with  comparison  with  measured  field  data  from  realistic 
testa.  The  numerical  procedure  to  test  the  general  validity  of 
this  approach  during  this  exploratory  phase  is  outlined  in  II-F 
and  la  Appendix  A.  Conclusion  of  the  time  allotted  for  the  numerical 
example  occurred  before  the  digital  computer  program  was  completely 
functional.  Complete  reliability  of  the  root- finder  subprograsi  to 
obtain  the  real  and  complex  roots  of  the  characteristic  equation 
was  not  achieved,  Two  methods  using  the  digital  program  and  another 
hand-calculated  test  method  were  used  in  searching  for  the  roots. 
Comparison  of  the  results  from  these  three  methods  lead  the  author 
to  conclude  that  a  combination  of  absolute  and  relative  root 
magnitudes  in  the  lower  modes  and  relative  root  magnitudes  in  the 
higher  modes  both  prevent  convergence  and  introduce  numerical 
Inaccuracies  into  the  root- finding  technique.  While  the  possibility 
of  programmer  procedural  error  is  not  precluded,  careful  checking 
has  revealed  no  blunders  (as  opposed  to  errors). 

The  major  program  using  a  double  precision  arithmetic 
method  for  automatic  computation  in  a  complex  number  mode  of  the 
pole  residues  and  tune  variance  appeared  to  function  properly. 
The  failure  to  achieve  general  program  reliability  naturally 
prevents  inclusion  of  the  velocity  responses  calculated.  However, 
for  the  modes  whose  characteristic  equation  roots  were  reasonably 
accurate,  responses  were  of  a  damped  oscillatory  form  as  would  be 
physically  expected.  As  the  weight  and  natural  frequency  of  the 
equipment  were  varied,  response  calculated  conformed  as  expected 
(25);  i.e..,  the  larger  the  mass,  the  lower  the  peak  mass  velocity 
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and  the  larger  the  rise  time.  The  influence  of  mass  of  equipment 
on  peak  velocity  can  he  seen  directly  in  equations  (2?)  and  (29). 
The  direction  of  the  mass  effect  on  the  rise  time  is  less  obvious, 
hut  it  is  clear  that  idle  mass  influence  is  directly  proportional  to 
its  magnitude.  The  accuracy  of  evaluation  of  the  incomplete  Bessel 
function  for  large  negative  root  values  during  early  times  of  the 
transit  period  became  progressively  worse  as  the  mode  order  increased. 
The  cause  of  this  13  shared  by  the  massless  shell  used  to  simplify 
root- finding  problems  and  the  form  of  the  equation  for  small 
times  involving  smell  differences  between  large  numbers*  Accurate 
results  for  these  early  phases  of  the  response  notion  can  be 
achieved  by  series  expansion  of  the  impedance  in  s  and  subsequent 
term  by  term  integration.  The  quasi-inpedance  would  then  be 
described  by  a  two-phase  function,  each  phase  with  a  time  validity 
range  The  need  for  this  added  complexity  needs  to  be  evaluated. 

Equation  (27)  describes  the  velocity  response  relative  to 
a  set  of  fixed  coordinates ,   wh  ch  is  adequate  for  seme  parametric 
studies.  However  ^  ic.  nest  situations  the  damage  evaluation  or 
"performance  level"  is  better  described  in  terms  of  response 
motions  relative  to  the  tall,  The  utility  of  the:  impedance  approach 
is  demonstrated  by  the  simplicity  with  which  the  relative  velocity 
is  obtained,  as  shown  by  equations  (32)  and  (33 ). 

In  general,  the  investigator  is  convinced  that  the  system 
model  is  ultimately  capable  of  providing  a  qualitative  and  quantitative 
estimate  of  the  resjjonse  of  a  real  cross  section  over  a  range  of 
test  conditions.  It  provides  an  analytic  model  bearing  some  f   b 


. 


■ 

J:".  •::.  . . 


..Lata  oe  * 


-33- 


charact eristics  of  a  distributed  mass  system  with  the  eonvenien 
of  a  discrete  Sjape&aoce  chain  system.  As  such  it  provides  a  means 
of  attaching  quantitative  meaning,  in  the  form  of  time  response,  to: 

(1)  shock- spectrum  nail  seen  at  equipment  fixed-base 
natural  frequency  in  discrete  mass  experiments 
investigating  distributed  mass  phenomena. 

(2)  knowledge  that  shock  spectrum  -  normal  mode 
approach ,  by  destroying  time  history  of  input  mat! 
results  in  over-conservative  estimated  response  (due 
to  neglecting  of  modal  time;  to  peak  in  summation 
process ) . 

Jty-  describing  the  response  in  terms  of  physically  meaningful 
parameters.,  the  author  believes  that  this  approach  will  allow  a 
better  understanding  of  the  basic  phenomena.  In  providing  a  method 
of  breaking  down  the  complicated  system  interaction  into  component 
interactions  capable  cf  description  by  engineering  approximation, 
the  impedance  approach  is  of  great  value. 
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1.  It  is  possible  to  formulate  an  analytic  approach  describing 
the  water  applied  pressure  loading  -  bull  -  equipment  inter- 
action with  the  use  of  engineering  assumptions,  capable  of 
solution  with  a  reasonable  amount  of  effort. 

2.  The  use  of  mechanical  impedance  and  transform  methods  provides 
a  powerful  analytic  tool  and  a  flexible  format,  aad  procedure 
for  the  introduction  of  more  exac^  model  components  and  of 
empirical  data. 
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rscqmmendations 

1.  Systematic  parameter  variation  ef  the  derived  quasi- impedance 
function  should  be  conducted  to  provide  rational  basis  for 
the  establishment  of  design  inputs  and  goals. 

2.  The  impedance  function  should  be  expanded  co  include  the 
effects  of: 

a.  multi-degree  of  freedom  equipment  -  foundation  model 

b.  foundation  and  equipment  damping 

c.  longitudinal  whipping  response 
d»  experimental  -  empirical  inputs 

3*  Response  variation  should  be  investigated  as  the  equipment 

attachment  point  varies  around  the  periphery. 
k.     Investigate  need  for  two-phase  impedance  function  to  extend 

accuracy  during  early  times  of  transit, 

5.  In  conjunction  with  2c  and  3,  investigate  the  effect  of  attack 
angle  of  incidence. 

6.  A  quantitative  investigation  of  the  'spectrum  dip'  phenomenon 
for  a  distributed  mass  model  of  the  form  developed  here 
should  be  conducted. 
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APPEHDIX 
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AR-PKNDIX  A. 
DETAILS  OF  PROCEDURE 

1.  Coa^puter  Procedur 

a.  In  order  be  expedite  tegi     la  t  ion  checking  of  -he 
IBM  7090  computer  program  daring  the  exploratory 
phase  of  the  problem,  the  order  of  the  characteristic 
equ&tioii  A  is  reduced  by  neglecting  the  shell  mas.*. 
As  a  result,  A  reduces  to: 

2  .  a  -2  m   .  a  -2 
oc2  a    oc2  * 

or 

2  _ 

»  "  a  s  +  ?a  g^   uhere 

rr  la  the  -loa-dljiienaional  shell  stiffness 
n 

-2  a 

n   *s 

pc 

Equation  (27)  becomes: 

i  „ a  ^  po  [(g + **k*m + *  K^y* 

^       (s>B)(a,,  8"  +  a_  s-*  +  a     s         a.,   3   '   a  ) 
*'xu  3  £  l  ° 

A2      ,.2 

titeHV      a^  *     c   •  g-  ♦  -75 


Ml      c2  c       1 


a     -  8«  .  ^  .  ~  ,     *  '  *1   '  a  *  ST 

"  a  J. 


Al      c2 

8 


'  Al   *  a   '  «q   '  q 


E- 

.1 

- 

fit  A  «A 


- 


•to 


-. - 


:as«*>d0tf  (YS)  aol&MtS 

■■}  a 


X* 


M 


—1*0- 


a.l  "  t*    n 


S.   -'T  -'C  0L.  &L, 


Rearranging,  with 

2     it 

jjp    natural  frequency  of 

2 


"J 


a2  2 

f    F   ,  —  -v 
°    °     pc3 


yields 

V  ~.  ... "■  ■■■■»  7  — — — a — -'-  -'%■■ »,    lA-l; 

"^   (a+B)  j  &  4  b,  a*5  +  b0  a*  +  b,  a  +  b 

L'      4|       *>       *■      °  : 

*2 


b2  =  [*»  «n  +  IJ  4ch  *n  *  ^  ykachJ 


-    a2  2 

Dl  Tj  ranch  cn 
c 


2 
bo  =  ^  <4eh  "n  «a 


The  root  .3  to  the  charge  '^eristic  equation  above  will 
be  call  x  ,   .   3 }   3  ,  3l< 

Inversion  of  (A-l)  ia  accomplished  by  calculating 
the  reaiouee  about  each  oi      -  \,   a,- 
Tbe  general  form  of  the  ti&f:  response  for  ail  modea 
other  than  translational3  i 
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iV      j. 
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«»ftOtt 
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v    (•)      v        t,\     ("  -a-B)V  (-B)  -  £  (l(a-X)     ('3>  >.Btt-D 


o  -a  v  3'  ,_(S  +  bjcs^  ■  3;.)(s'  -  80)(s-  -  s,.)_r 
vo  Hi  <V  Wpy  Wb^;  ^R«^  -  *£hsJ  *  »3)-T 

(a  \  i/j  (A-2) 

b.  Shell  and  equipment  -  foundation  parameters  were 
introduced  as  data  to  the  machine,  which  calculated 
the  modal  coefficients  of  the  characteristic  equation 
of  (A-l)  for  the  first  twenty  modes,  Subroutine 
root-finder  was  programmed  to  search  for  the  real- 
complex  roots  of  the  resulting  qu&rtic.  The  first  sub- 
program used  was  aii  XBH  9MS8  Distribution  Program 
.-,  :"MJLU5B"*.  Persisting  machine  overflows,  I 
perhaps  to  coefficient  incompatibility,  resulted  in 
the  author  prograiaalag  the  Euler  -  Descartes  Met! 


*"MQLlM,r"is  based  upon  the  theory  of  Miller,  J.,  *A  Method  for  the 
Solution  of  Equations  with  Automatic  0         <ors",  Mathematical  Tables 
of  Autcg^tlcComjtatatioa;  1956,  pp  208,  &f« 
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for  closed  form  solution  of  quartic,,  previously 
used  to  check  computer  results.  The  basic  steps  are; 
1.  reduce  quartic 

k  3    2 

y  +  pyJ  +  qy  +  ry  +  s  =  0 

to  form 

k  2 

x  +  ax  -*  hx  -t-  c  =  0 

by  substitution 

y  •  (X  -  p/4). 

We  find 

2 


a  ,  6* :  ft' 


,   l6r  -  8pq  -;-  2p* 

-3p4  +  I6gp2  -  64rp  +  2^6  s 

c         »^ 


2.  With  i,  m,  and  n  as  the  roots  of  the 
resolvent  cubic; 

+3  J    /»\  *2  „  /&2  -  ^\  +   0    r» 
tJ  +  C^)  t  +  C — jg ;  t  -  ^  «=  0 

The  reduced  quartic  roots  are: 


;;. 


x2  •  ^  (-  fl     ->-  fm      ♦  /~a~    ) 

where  upper  sign  is  used  if  b  is  positive 
and  the  lover  sign  if  b  is  negative. 
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c.  Having  calculated  the  roots  of  the  characteristic 
equation,  the  computer  was  programmed  to  evaluate 
the  modal  incomplete  Bessel  function  I.  of  equation 
(B-l)  at  three-degree  increments  in  envelopment  angle 
with  a  three-staticu  Simpson's  Rule*  After  storing 
these  values,  the  residues  of  equation  (A-2)  were 
evaluated  at  intervals  of  the  non-dimensional  time  T 
corresponding  to  the  envelopment  angle  (0)  increments 
of  the  integrated  Bessel  function* 

2.  Typical  Shell  and  Equipment  -  Foundation  Parameters  Used. 
a.  Using  typical  values  for  E,  I,  a,  p,  p  ,  m  , 

3    IX 

yielded  for  the  simple  shell, 

<a  *  (extensions!.)  -  (1.13  x  10°)  £=~, 

sec 

i *   (flexuraf )  -  **  ^  '/f  x  782  Sfl*  , 

1  +  n         sec 

^  2  »  4.23  x  10      lh../ft2  surface/ ft  length 
o  I 

<\2  ■  3035.   (1  -  n2)2    lhf/ft2  surface/ft  length. 

fc.  Three  combinations  were  chosen  for  the  equipment  - 
foundation  parameters: 


1.     K^  «  155.2  slugs/ft. 

k    =  5.515  x  106  lbf/ft. 

f    «=  30  cycles  per  second  (cps) 


2.     K^  =  7-76  slugs/ft 
k    ■  I.960  x 
f      =80  cps 


k    -  I.96O  x  106  lbf/ft. 


■ 

T    6. 

.a 

- 


• 
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f  -  12-0  cps 

c>  Saoek  Wave  Characteristics. 

A  peak  prejsur«  of  2000  psi  aad  a  decay  constant 
0  »  20.8  see"  were  introduced  as  data  to  the 
machine . 


. 
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UUMERICAL  EVALUATION  OF  IHCGMi'LEIE 
AND  COMPLETE  MODIFIED  W3SEL   F11HCTION 
OF  THE  PIR3T  KIHT 


During  the  shock  wave  transit  across  the  cylinder  <, 
evaluation  of  the  Fourier  moo,  i  coefficients  of  the  incident  pressure 
yields  the  function: 

In,  (s)  -  i  '  eB  cos  8  cos  no  d  (b-1) 

o 

0  «   «  2n 

where  0  is  the  angle  of  envelopment.  For  times  greater  than  the 

transit  time  of  the  shock  wave.  I  m  becomes: 

n9 

an 

Ia  (•)  ■  |  J  e8  cos  °  cos  n9  do  .  (B-2) 

b 

Equation  (C-2)  is  the  definition  of  the  modified  Bessel  function 

of  order  n.  Equation  (c-1)  is  therefore  termed  the  incomplete 

modified  Bessel  function  of  order  n.  Equation  (C-l)  is  not  integrable 

in  closed  form^  but  is  capable  of  evaluation  by: 

(1)  numerical  quadrature 

(2)  rapidly  convergent  series  approximation. 

Figure  B-l  plots  I  .  (s)  as  function  of  T  (or  equivalently  )  for 
a  specified  real  root  of  the  characteristic  equation,  as  obtained 
using  three- station  Simpson's  Rule,  It  is  presented  to  give  an 
idea  of  the  gross  form  of  the  function.  In  general,  s  is  complex 
and  requires  integration  of  real  and  imaginary  parts.  Equation  (B-l) 
takes  the  form: 


(X-  *o:>  V;     £  «  (8) 

3918  aaatf*  io*t      .JtastigaXaTO  < -gee  axtt  e.;  8  yiadv 

: aaawwrf      _  .o  «al^  ttanai* 

(S-«)  .  a  «oo  ■'  •^  ••  ]     i  .  (a)  JL 

aattasrt  Isaaaff  feJSlftos?  ad*  to  acJttlails*  satt  ai  (S-o)  aolJaayS 
^»oai  sm*J  JbwKKte,  .1.-^  ei  *xpl     .a  tmirto  to 

»£da««a  •!  Jwf  <snol  toaold  at 

*ixtf*l*jw/p  l*-.aiiaia.xi     (I) 

xtxaaqi**  e  jiisfi  jco^isvxioa 

UmSjn  .  iomA  aa  (a)      i  ».  ~f  aixi&n 

**  *£  .inrfaaiaila  aaH  i  Lbs?  aallxoaqa  a 

xte  ©vis  oS  betMatQKq  ax  *I     .six*  a 'oca.  -aewtt  s&taxr 

Xanuag  m     .a»i  .-.uoig   >xu-  ip  Mbx 

iI*€'  nc  •••  -Jmlimt  jorui  J>  .ol#,ir%>«at  aaiJUf>ai  Aob 

:«tot  ajbr  asXrt 


e 


Ine  (s)  =  i     f  e(x+ly)  cos      coa  ne  |«  ,  which  can  toe 


0 
written  by  Baler's  identity  as: 

9 


j  e  3      ,  cos  (y  cos  )  *  1  sin  (y  cos  )  cos  n  I 
0 


1 


X  cos 


e      cos(y  cos  )  cos  n  il 


0 


+  1 1  |  ex  coa  ?   sin  (y  cos  r )  cos  al  II        (B-3) 
0 

waere  x  and  y  refer  to  the  real  cad  imaginary  parts  of  s  respectively. 

A  digital  machine  mkes   the  evaluation  of  (B-3)  >ite  effortless. 

An  alternate  method  of  evaluation  is  to  use  «he  Jacobi-Anger 

formula  and  a  series  expansion  of  the  modified  Bessel  function. 

Bes3el  functions  of  non-negative  order  n  =  0,  1,  2. . . .  are  single 

valued  integral  functions  of  s.  They  may  toe  formed  as  coefficients 

of  the  Fourier  series 

n°*» 

c+  z  cos  %m  \     ^  ^  ^  co3  a  t      (jacotoi- Anger  Formula) 

n=0  (B-4) 

where  ^  (Neuaaaa  factor)  =  1  I tor  n  ■  0 

=  2  for  li  m  1,  2, . ••     . 
Using  (l-^)j    (B-l)  goes  to: 

n  Ine  (s)  *  icos  a      ^m  3^  (s)  cos  sir  d?  . 
0  i=o 


atf 


■ 


:  n  ■■    , 


.-?  m$*i 


KM  ( 


40  '  ~  I  * 


■ 


■r 


JMlV)         . 


-   ,  . 


total 


j 


0  o: 
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The  e     I  (a)  is  not  a  function  of  the  integration,  so  that  we  can 
bring  it  outside  the  summation  sign  and  interchange  the  integration 
and  summation  signs  yielding  a  form: 

0  e 

rr  I  (s)  ld(s)  J  cos  no  It  +  3  1^  (s)  j  cos  M  cos  I  48 

0  0 

e 
+  2I2  (s)  j  cos  n*  cos  29  de 
0 

e 

*  ft*  (s)   cos  n  cos  3  d 
0 

+♦..♦.   •  (B-p) 

Equation  (B-5)  involves  for  real  s  only  tabulated  functions  and 

integrations  soluable  in  closed  form>  Further,  because  of  rapid 

relative  decay  of  the  modified  function  with  increasing  order, 

five  or  six  terms  should  be  adequate  for  engineering  accuracy. 

For  non-tabulated  arguments  the  modified  Bessel 
function  of  the  first  kind, 

I  (s)  •>  i"m  J  (is), 

where  J  (a)  is  the  Bessel  function  of  the  first  kind  is  defined 
m 


by: 


.  M.J    IxZafUJL  (B-6) 

m  2  h 

S /-,      3      ^        S 

'  2m  r  ,  +1v  KX       2(2m+2)   2.M2m+2)l2a+^) 

+  ...)     (  l&BG  s|  <  tt) 
with  r  (n)   gamma  function  ■  (n-l)!  ■  (n-l)  V   (n-l). 


(a)     I 


tfl  ~'v..:.     | 


S  * 


0 


- 


.( .     I 


.<;    . .  ■;:  -v     ..lis    :.o    r.\..  i 

- 


<(«- 


■ 


W    ** 


■ 


-' „» 


-48- 


For  •mall  values  of  s,  (B-6)  can  be  approximated  by  the  leading 
terms: 

vw'A':'  (B-7) 

For  large  values  of  9, 

im  <•)  =  ~^r;  .  (B-8) 

In  manipulating  the  modified  Beasel  functions^  the  following  re- 
currence relations  are  useful, 

h  la  (a  8)  =  a  Vi (a  8)  -  f  ^  (a  8)  ^^ 

h  Xa  (a  8)  *  a  Vl  <a  3>  *  !  Xa  <a  ">  <B-10> 

2  S  ^.  <a  •)  -  a  [V*  (a  S)  +  Vl  {a  S)        (B"al) 

Vi  <a  s>  -  Vi (a  8)  -  H  ^ (a  b)  •  (B-12) 
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apfehpix  c. 

FOUR-POLE  PARAMETERS 


1.  Development  of  Four-Pole  Parameters  for  Fundamental  Elements. 


a. 

Mass 

*i* 

MUM 

F   •  V 

1 

if 

Perform 

ance  Equations  arc: 

v"    e    V       m  t 

vl       v2 

F     -    F 
*1        *2 

=    8Illf 

where  F  and  V  are  arbitrary  complex  functions  of 

time  and  F  and  v1  are  the  transforms  of  these  functions. 

It  is  clear  that  tho  four-pole  matrix  is  then: 

fal  -  [*  J*] 
h.  Massless  Spring  Four  Poles 


«i«  v5 


3PRIH0 


*2'  V2 


2 


Performance  Equations  are: 


F  -  Fx  -   P2 


or  V  =-?  +  V 


The  four-^ole  matrix  is  then 


£*i  .- 
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s/k 
Viscous  Damper 
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vvi 


DAMPER 
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\mm  mc 

i     . 

*     •    % 

•  XV)  j  -     *]  *  «  xi 


•xo 
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FerftattHM  txons  are: 

c(vx  *  v2)  or 


1       C     2 


r2   ' 


The  four-pole  matrix  is  tc. 

m  •  &  a 

2-  Combinations  of  Component  Quadripoles; 
a.  Series  Combination. 

If  the  output  of  one  component  tm   sonnteted  to  the 
input  station  of  another  system,  the  two  systems  are 
in  series.  For  example,  if  the  two  components  ta 
are 


connected  at  point  2,  the  components  are  connected 
in  series.  If  the  four-pole  of  flj  is 

Lc  dJ 

the  four-pole  of  [2]  is 

four-pole  relating  F-  and  V^  to  F-  and  V    is 
■•"A    B"l  rs    m 

Lc  dJ  Li  e 

So  that  the  net  series  connected  system  four-pole 
is  equal  to  the  matrix  product  of  the  four-pole 
matrices  associated  with  each  of  the  systems  making 
up  the  composite  system. 


. 
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. 

fTut 

V* 

1 

3TB  ec 
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Ml 

■ 1 

si: 


»i  ^V  b*t  i^     ,  ^ 
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b.  Parallel  Combination. 

Components  are  considered  bo  be  connected  in  parallel 

(1)  all  the  input  and  output  Junctions  move 
with  the  same  velocity 

(2)  tiie  net  input/output  forae  of  the  resolvent 
quadripole  is  equal  to  the  sum  of  the  input/ 
output  forces  of  the  constituent  four  poles, 

If  the  system  obeys  Maxwell's  Reciprocity  Principle, 


vwmm 

V 

V    1 

"    V2 

w 

1 

F  ,  V 
2*     2 

2 

can  be  replaced  by 

h 

1 

1 

2 

llxe  asutrix  s  of  the  composite  system  relates  input 
to  output  by 

!_vr!   LKj  l\'2 


where  3 


I  »/b  -i/b"i 

'  Ll/B  -a/B.I 

JD^  -1^ 
ll/B,   -A,/Bj 


1?B,   -A0/B 


2^2  ' 


i 


■•   *u*jm£    NU 


-^^. 


.    of  t&e  abc  •  •  gpwilftg  to  a  apring- damper 
p&rali  .1  uorabiatation  yield  an  eqaivaleat  four-pole 

\X 

Irs+k        | . 

More  extensive  developments  are  presented  in  (23)> 
(1),  MM   {2k). 
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£CCJK 
IK-VACUO  NOB  I R  -tffmBW  SHELL 

Analytic  methods  are  available  providing  an  additional 
aaasure  of  realism  in  the  description  of  the  modes  aad  frequencies 
of  in-vacuo  riag-s&il'icaed  cylindrical  shell.     Because  of  the  ainrple 

extension  required  to  apply  these  frequencies  to  she  quasi-itapedanee 

(45) 
developed^  the  general  form  of  an  appro:  aethod  hy  Bleiehx  " 

tki) 

as  reported  in  v  J/  are  included  here,  Results  of  the  approximate 

tk  ) 
engineering  analysis  of  v  *//  yield  results  in  agreement  with  model 

experimenta  within  ♦  10  percent. 

IpMCBClature 

=  frequency  of  shell 

a).    «  frequency  of  ring 

a     =  mass  of  shell  per  unit  of  area 
s 

m  m  mags  of  ring  frame  per  unit  length 

L_     frame  spacing 

I     n  second  moment  of  inertia  of  ring  frame 

section  plus  a  length  of  shell  plating 

equal  to  2/p. 
K  &  K  *  Bleich  constants  tabulated  in  (46) 
o     *  thic     of  shell 
t>       oisson's  ratio 
E     *  Young's  Modulus 
a     ■  mean  shell  radius 
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Bleich  represents  $Mi  total  frequency  of  the  in- vacuo 
vibration  of  stiffened  cylinders  as  a  suxaaation  of  the  shell 


frequencies  plus  the  ring  frequencies. 
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where  m  «*  "_  x  1^,  x  a  . 
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